Tunable bandgap and magnetic ordering by adsorption of molecules on graphene 
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We have studied the electronic and magnetic properties of graphene and their modification due 
to the adsorption of water and other gas molecules. Water and gas molecules adsorbed on nanoscale 
graphene were found to play the role of defects which facilitate the tunability of the bandgap and 
allow us to control the magnetic ordering of localized states at the edges. The adsorbed molecules 
push the wavefunctions corresponding to a-spin (up) and /3-spin (down) states of graphene to the 
opposite (zigzag) edges. This breaks the sublattice and molecular point group symmetry that results 
in opening of a large bandgap. The efficiency of the wavefunction displacement depends strongly 
on the type of molecules adsorbed. 



Monolayer graphene [l| with its high electron mobility, 
unique magnetic phenomena [2j] , and unusual relativistic- 
like properties [3[ has generated an upsurge of activities 
in materials science. Materials exhibiting magnetic prop- 
erties are in great demand for applications in nanoscale 
electronics and spintronics. Most magnetic materials are 
metals, where ferromagnetism (which results from an 
imbalance between the spin-up and spin-down unpaired 
electrons) is often destroyed by thermal fluctuations. In 
carbon systems the magnetic properties are stable even 
at room temperatures. The nonmetalicity of carbon sys- 
tems makes them biocompatible and therefore are ideal 
for a wide range of applications not only in nanoelec- 
tronics and spintronics 0, H, @], but also in medicine 
and biology. In most carbon systems the origin of mag- 
netism is unclear, but in graphene it is expected to be 
the result of spin ordered states localized at the edges 



The absence of a bandgap is however a major hin- 
drance for graphene's application in nanoelectronics. 
Various mechanisms for opening a gap have been de- 
veloped that involve breaking of certain symmetries in 
graphene by defects , an a ppl ied bias [% [H[ , and in- 
teraction with other materials 11411 . Noting that the con- 
ductivity of the carbon systems [15|, 16] is extremely sen- 
sitive to adsorption of gas molecules due to the charge 
exchange between them, we assume that a more control- 
lable way to induce interactions that break a symmetry 
of graphene would be gas adsorption. The charge ex- 
change between the adsorbed molecules and graphite was 
found to be rather low [17[ , and donor or acceptor behav- 
ior is exhibited depending on the type of molecules ad- 
sorbed. However, since gas adsorption changes the elec- 
tronic properties of graphene due to the charge exchange 
between them, we suspect that adsorption should affect 
the localization of electronic states along the edges and 
can facilitate opening of a gap. 

Here we study the influence of adsorption of water on 
the electronic and magnetic properties of graphene by 
molecular mechanical methods using the spin-polarized 
density functional theory with semilocal gradient cor- 
rected functional (UB3LYP/6-31G 0) in the Jaguar 
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FIG. 1: (color online). Nanoscale graphene structures, where 
N is the number of rings along the zigzag edges and m is the 
along the armchair edges. 



6.5 program [19[. The van der Waals interactions, which 
impact the interaction between graphene and adsorpant, 
are not considered in DFT, that lead to underestimation 
of the adsorption energy [l7[. However, the DFT is reli- 
able for investigation of the influence of the adsorption on 
the electronic properties of the graphene. We chose the 
nanoscale graphene structures as shown in Fig. 1. The 
number of rings along the armchair edges was varied as 
m = ... k or m = 1 ... (k — 1) to take into account 
different geometries of the graphene seams (indicated in 
Fig. 1), where k was always odd. The number of carbon 
rings along the zigzag edges N was varied from 3 to 7. 

For a nanoscale graphene the highest possible symme- 
try is D2h, i.e., the planar symmetry with an inversion 
center, where all four seams are characterized by the same 
boundary conditions and exhibit the lowest energy in the 
system due to the confinement effect. The optimization 
process makes these seams structurally identical with the 
same density distribution. Therefore, the seams have the 
same spin order and play the role of the center for spin or- 
dering along the edges (see Fig. 2(a)). The resulting spin 
polarization exhibits a ferromagnetic ordering between 
the seams and between the zigzag edges, but a mixed or- 
dering along the zigzag edges. Most molecular orbitals 
including the highest occupied and the lowest unoccupied 
(HOMO and LUMO) exhibit zero spin density, because 
they have the same density distribution for the a- and j3- 
spin states, but with opposite sign. For the HOMO and 
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LUMO each spin state is equally delocalized over the two 
zigzag edges. The molecular orbitals contributing to a 
non-zero spin density are located far from HOMO-LUMO 
gap and they are not localized along zigzag edges. The 
spin density and spin polarization for graphene of size 
iV=4 are presented in Fig. 2 (a). For nanoscale graphene 
optimized without any symmetry, with the planar sym- 
metry CI or with the C2v symmetry, the seams are not 
constrained to be identical. That makes the zigzag edges 
exhibit the lowest energy and as a result, a ferromag- 
netic ordering of the spins along the zigzag edges with 
a bandgap of ~ 1.6 eV for C2v symmetry and CI sym- 
metry, instead of ~ 0.5 eV for D2h symmetry. In these 
cases, the HOMO and LUMO are characterized by the a- 
and /3-spin states localized on the opposite zigzag edges. 
Although the state of graphene with D2h symmetry is 
comparable in energy to the states of C2v and CI sym- 
metries, the difference between them is already less then 
0.02 eV for the N >6 and m > 7. This implements the 
high metastability of the ground state of C2v symmetry. 
Therefore, D2h is chosen here in order to study the in- 
fluence of adsorption of gas and water molecules on both 
spin ordering and the bandgap. 

The D2h symmetric graphene structures with an even 
number of rings along the zigzag edges N = 4, 6 and the 
number of rings at the armchair edges that are varied as 
m = . . . fc, are found to be the only case with a non- 
zero spin density (Fig. 2 (a)). For these structures and 
sizes, the atomic wave functions from the opposite seams 
cancel each other at the center of the edge, where they are 
out of phase. In that case, the a- and /3-spin states are 
spatially separated which leads to non-zero spin density 
for certain orbitals. The a state is localized on the seams, 
while the /3 state at the center of the zigzag edges (see 
Figure 2 (a)). We found this effect for several molecular 
orbitals that are far from the HOMO-LUMO gap, while 
other orbitals including the HOMO and LUMO have the 
zero spin density. These states are non-degenerate, that 
diminish with increasing structure size. For example, the 
energy difference between the a- and /3-spin states for 
HOMO is 0.0764 eV if 7V=4 and 0.00625 eV if N=6. This 
can explain the experimental observation of the single 
local density peak in the valence or the conduction band 
for large graphene structures 20|, Hi]. For an odd number 
of carbon rings along the zigzag edges N = 3, 5, 7 and 
m = . . . fc, the spin density at the edges is close to zero 
because the a- and (3- degenerate states exhibit almost 
the same electron density distribution over all molecular 
orbitals. In a structure where the armchair edge is varied 
as m = 1 ... (fc — 1), we found the opposite behavior of 
the spin degeneracy: Degenerate a and (3 spin states are 
observed for the even ring number on the zigzag edges, 
while for odd N the states are non-degenerate. 

Our most striking results are for the case of adsorbed 
water molecules on nanoscale graphene. The interaction 
between water and graphene is weak because the attrac- 



tive forces between the water molecules within a cluster 
are stronger than the attractive force of the graphene 
surface. Therefore, several water molecules will tend to 
form a water cluster, that was studied for graphene struc- 
ture of sizes m = . . . k and 7V=3,4,5,6. A single water 
molecule on the graphene surface is localized at the cen- 
ter of a carbon ring, ~ 3.0 A above the graphene sur- 
face. The optimized position of the water molecule on 
the graphene surface corresponds to its dipole moment 
being perpendicular to the graphene surface and directed 
inward, which provides an effective charge donation from 
the oxygen to the graphene through the unshared pair 
of p electrons of the oxygen. When we increase the 
number of water molecules, this water molecule becomes 
a "cluster link" to the graphene, i.e., the extra water 
molecules are connected to the cluster link by the hydro- 
gen bonds. Therefore, the charge donation to graphene 
occurs mostly through oxygen in the cluster link, while 
acceptance by the cluster is through the hydrogen atoms. 
For the graphene structure with N < 4, where the con- 
finement effect is strong, the interaction between the arm- 
chair edge and the water cluster leads to unrolling of the 
water cluster in such a way that one termination side is 
attached to the armchair edge and the other - the cluster 
link - to the center of a carbon hexagon. The attachment 
to the armchair edge is a result of its low energy in com- 
parison to the zigzag edge [22]. The reconstruction of the 
zigzag edge can lead to lowering of its energy that would 
change the direction of cluster unrollment. For N > 4 
the water cluster is formed at the center of graphene. 

Water can be treated as a defect on the otherwise 
defect-free graphene surface. Figure 2 (b) shows that this 
defect repels the wavefunctions of the molecular orbitals 
corresponding to a- and /3-spin states towards opposite 
zigzag edges of the nanoscale graphene. This leads to spin 
ordering of the localized electrons on the edges such that 
the a-spin state is localized on one sublattice while the 
/3-spin state on the other sublattice. The non-zero spin 
density and repulsion of the wavefunctions by water were 
observed for all types of graphene structures independent 
of size. Just as for pure graphene, the size is found to be 
crucial only to preserve the non-degeneracy of the a- and 
/3-spin states. The defects affect the wavefunction distri- 
bution only for the nearest carbon rings. This means 
that, depending on the graphene size and the size of the 
water cluster, the wavefunctions can be repelled locally, 
or over the whole structure. For a graphene size of N = 6, 
the water molecules build a cluster close to the graphene 
center. Then the a- and /3-spin wavefunctions are pushed 
toward the edges only for the localized states along the 
line of the water location (Fig. 3). 

The repulsion of the a- and /3-spin wavefunctions to op- 
posite zigzag edges discussed above, breaks the molecular 
and sublattice symmetries and opens a gap. The magni- 
tude of the gap directly depends on the interaction of the 
graphene surface with water, particularly on the charge 
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FIG. 2: (color online). Spin polarizations in nanoscale graphene of size iV=4 and m = . . . k: The spin density and electron 
density for the a-spin (up) and /3-spin (down) states are plotted for isovalues of ±0.01 e/A 3 . For the spin density, different 
colors indicate the a- and /3-spin states, while for the electron density distribution (a- and (3- states) the different colors 
corresponds to different signs of the molecular orbital lobes, (a) The electron density distribution in pure graphene for a-spin 
(^homo-4 = —7.119 eV) and /3-spin (Ehomo-4 = —7.295 eV) states for the fourth orbital below the highest occupied molecular 
orbital (HOMO-4). In pure graphene, a non-zero spin density is obtained only for an even number of carbon rings along the 
zigzag edge (N = 4,6, . . .). (b) The electron density distribution of the a-spin state (Ehomo = —5.386 eV) and the /3-spin 
state (Ehomo = —5.270 eV) for the HOMO for graphene with adsorbed water. The water cluster plays the role of a defect 
pushing the a- and /?-spin wavefunctions to the opposite edges. 
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FIG. 3: (color online). Spin polarizations in nanoscale 
graphene of size N=6 and m — . . . k. The spin densities 
are plotted for isovalues of ±0.01 e/A 3 : (a) in pure graphene; 
(b) for graphene with adsorbed water. 



exchange and the association energy. For a single water 
molecule, if the dipole moment of the water molecule is 
directed inward (toward the graphene surface), the gap 
increases by up to ~ 2.0 eV, but if the dipole moment 
is directed outward then the charge donation from oxy- 
gen to graphene does not occur and the gap opens by 
only up to 0.8 eV. For the inward directed dipole, the 
gap slightly decreases with increasing distance between 
the water molecule and the graphene surface. For the 
water cluster, the cluster link usually is a donor and its 
interaction with the graphene surface defines the associa- 
tion of the water cluster to graphene and the efficiency of 
opening a gap. The orientation of the cluster link to the 
graphene surface is changed due to the intrinsic interac- 
tion within the water cluster. Therefore, the association 
energy of the water cluster on graphene and hence the 



graphene gap oscillates depending on the orientation of 
the cluster link to the graphene surface and correlation 
between attractive forces inside the water cluster and at- 
tractive forces of the graphene surface. For a cluster con- 
taining three water molecules, the association energy has 
a global minimum and the same minimum is observed 
for the graphene gap (Fig. 4 (a)). With an increase of 
the graphene size, the gap decreases (~ 1.2 eV for N > 
6) due to the vanishing of the confinement effect and 
the gap oscillation diminishes (the dashed lines in Fig. 4 
(a)). The energy gap of 1.2 eV is in perfect agreement 
with earlier investigation of the bandgap of nanoribbons 
with ferromagnetic ordering along the zigzag edges by 
the B3LYP method [5|. We predict that using a water 
detergent would be a way to reduce the attractive forces 
inside the water cluster and thereby increase the water- 
graphene interactions. This would allow one to tune the 
energy gap in graphene and improve the repulsion of the 
wavefunctions towards the opposite zigzag edges. 

Opening of the gap is a result of breaking the sublattice 
and molecular symmetries of graphene, but the gap size 
can vary depending on the interaction between the adsor- 
pant and the graphene surface. Therefore, we also inves- 
tigated the possibility to tune the energy gap by changing 
the type of adsorbed molecules. The HOMO-LUMO gap 
for nanoscale graphene with adsorbed NH3, H2O, HF, 
or CO molecules are presented in Fig. 4 (b). Since the 
gap opens due to the symmetry breaking, the interac- 
tion of nanoscale graphene with any molecule leads to a 
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FIG. 4: (color online). The effect of adsorption on the gap 
in nanoscale graphene. (a) The HOMO and LUMO energies 
for graphene of size N=3 (solid line) and N=5 (dashed line) 
versus the number of adsorbed water molecules, (b) Influence 
of adsorption of different adsorbed molecules on the bandgap 
of the nanoscale graphene of size N=4. 



gap opening for up to ~ 1.5-2.0 eV. The ability of the 
adsorbed molecules to push the wavefunctions of a- and 
/3-spin states to edges is found to differ between the type 
of molecules chosen and the HF and CO were least effec- 
tive. Moreover, the HF adsorption leads to spin-selective 
semiconductor behavior of the graphene. 

In summary, the adsorbed molecules on graphene are 
found to act as defect centers which alters the magnetism 
through the repulsion of the wavefunctions correspond- 
ing to a- and /3-spin states to the opposite zigzag edges. 
This results in a bandgap opening and stabilization of the 
state with parallel spin ordering along the zigzag edges 
and anti-parallel one between the opposite edges. The 
reversible breaking of the graphene symmetry and open- 
ing of a gap due to adsorption, as obsevred in experi- 
ments 
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, would open a new direction in control- 
ling the electronic properties of graphene. In particu- 
lar, transition from metallic to semiconductor state by 
molecular adsorption is a promising route to develop fast 
on-off switching devices [23|, super-sensitive gas sensors 
[Til . [l6[ and FET transistor [24]], where adsorption can 
help ambipolar-to-unipolar conversion. For FET appli- 
cation, the graphene can operate similar to carbon nan- 
otubes. The gap opening in graphene induced by adsorp- 
tion can provide ambipolar- unipolar conversion of FET, 
as a result of occurrence of a significant Schottky bar- 
rier at metal/graphene contact, leading to modification 
of the conducting channel operation. The spin-polarized 
state of graphene stabilized by adsorption provides an 
opportunity to manipulate spin current by an in-plane 
electric field applied across the zigzag edges Q • The ap- 
plied field increases a bandgap for one-spin state, while 
suppressing the other, thereby generating half-metallicity 
in graphene. This effect can be used for developing spin- 
related electronic devices 0, 0, S , such as a spin- valve 
[6] and digital memory devices |25j . 



* Electronic address: chakrabt@cc.umanitoba.ca 
[1] A.K. Geim and K.S. Novoselov, Nature Materials 6, 183 
(2007). 

[2] Y. Zhang, Y.-W. Tan, H.L. Stormer, and P.Kim, Nature 

438, 201 (2005). 
[3] V.M. Apalkov and T. Chakraborty, Phys. Rev. Lett. 
97, 126801 (2006); X.-F. Wang and T. Chakraborty, 
Phys. Rev. B 75, 041404 (2007); M.I. Katsnelson, K.S. 
Novoselov, and A.K. Geim, Nature Phys.2, 620 (2006); 
H.-Y. Chen, V.M. Apalkov and T. Chakraborty, Phys. 
Rev. Lett. 98, 186803 (2007). 

Y.-W. Son, M. L. Cohen, and S.G. Louie, Nature 444, 
347 (2006). 

E. Rudberg, P. Salek, and Y. Luo, Nano Letters 7, 2211 
(2007). 

S. Cho, Y. Chen, Y. and M.S. Fuhrer, Appl. Phys. Lett. 
91, 123105 (2007) 

Y.-W. Son,M. L. Cohen, and S.G. Louie, Phys. Rev. Lett. 
97, 216803 (2006). 

M.Y. Han, B. Ozyilmaz, Y. Zhang, P. Kim, Phys. Rev. 
Lett 98, 206805 (2007). 

K. Nakada, et. al., Phys. Rev. B 54, 17954 (1996). 
H. Lee, Y.-W. Son, N. Park, S. Han, and J. Yu, Phys. 
Rev. B 72, 174431 (2006). 

L. Pisani, B. Montanari, and N.M. Harrison, New Jour- 
nal of Physics 10, 033002 (2008). 

F. Cervantes- Sodi, G. Csanyi, S. Piscanec, A.C. Ferrari, 
Phys. Rev. B 77, 165427 (2008). 

M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusak- 
abe, J. Phys. Soc. J. 65, 19203 (1996). 
S.Y. Zhou, D.A. Siegel, A.V. Fedorov, A. Lanzara, Phys. 
Rev. Lett. 101, 086402 (2008). 

F. Schedin, A.K. Geim, S.V. Morozov, E.W. Hill, P. 
Blake, M.I. Katsnelson and K.S. Novoselov, Nature Ma- 
terials 6, 652 (2007). 

T.O. Wehling, K.S. Novoselov, S.V. Morozov, E.E. 
Vdovin, M.L Katsnelson, A.K. Geim, and A.I. Lichten- 
stein, Nano Letters 8, 173 (2008). 

F. Ortmann, W.G. Schmidt and F. Bechstedt, Phys. Rev. 
Lett. 95, 186101 (2005). 

Becke A. D., J. Chem. Phys.98, 5648 (1993). 
Jaguar, version 6.5. Schrodinger. LLC: New York, NY, 
2005. 

Z. Klusek, Z. Waqar, E. Denisov, T. Kompaniets, I. 
Makarenko, A. Titkov, and A. Bhatti, Appl. Surf. Sci. 
161, 508 (2000). 

Y. Kobayashi, K. Fukui, T. Enoki, K. Kusakabe, and Y. 
Kaburagi, Phys. Rev. B. 71, 193406 (2005). 
P. Koskinen, S. Malola and H. Hakkinen, Phys. Rev. Lett 
101, 115502 (2008). 

X. Dong, et. al., Phys. Rev. Lett. 102, 135501 (2009). 
Y.-M. Lin, J. Appenzeller and P. Avouris, Nano Lett. 4, 
947 (2004); D. Kang, et. al., Appl. Phys. Lett. 86, 093105 
(2005). 

[25] D. Gunlycke, et. al., Nano Lett. 7, 3608 (2007). 



[4] 
[5] 
[6] 

[8 

[9; 
[10 

[11 
[12 
[13 
[14 
[15 

[16 

[17 

[18 
[19 

[20; 

[21 

[22; 

[23 
[24 



